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^ '. Abstract 

O 

r-| i We consider the dynamical system consisting of a quantum degree of freedom A 

interacting with A^ quantum oscillators described by the Lagrangian 



^: r- ii2^^/1^2_l, 



L = lA^ + J:[lx!-lim^ + e^A^)x^] 



In the limit A^ -^ oo, with e'^N fixed, the quantum fluctuations in A are of order 1/A^. 
In this limit, the x oscillators behave as harmonic oscillators with a time dependent 
mass determined by the solution of a semiclassical equation for the expectation value 
{A{t)). This system can be described, when {x{t)) = 0, by a classical Hamiltonian 
for the variables G(t) = (x^(t)), G(t), Ac(t) = {A(t)), and Ac(t). The dynamics of 
this latter system turns out to be chaotic. We propose to study the nature of this 
large- A^ limit by considering both the exact quantum system as well as by studying 
an expansion in powers of 1/A^ for the equations of motion using the closed time path 
formalism of quantum dynamics. 



1 Introduction 

The definition and observation of cliaotic beliavior in classical systems is familiar and 
more or less well understood [IH . However the proper definition of chaos for quantum 
systems and its experimental manifestations are still unclear 0. Here we first study 
a simple model of two coupled systems which displays semiquantum chaos when 
one of the systems can be treated "semiclassically." We then study a purely quantum 
system of A^ + 1 degrees of fredom which has the identical dynamics in the large 
N limit as our original system. In this way we can determine, as a function of A^, 
what is the time scale for quantum fluctuations of the "classical" oscillator to be of 
significant size. We also can determine how this time scale is related to the time scale 
determined by the maximum Lyapunov index. The question we are interested in here 
is whether quantum fluctuations become significant before or after the original system 
is sensitive to initial conditions. 

In a classical chaotic system, such as the weather, we are accustomed to situations 
where there is lack of long time forecasting because of the sensitivity of the system 
to initial conditions. The simple model system considered here has the unusual fea- 
ture that one has to give up long term forecasting even for the quantum mechanical 
probabilities, as exemplified by the average number of quanta at later times [^ . 

First let us review the original system which displayed semiquantum chaos. Con- 
sider two coupled quantum systems described by the Lagrangian, 

L = h^ + \A'-\{m^ + e'A^)x^ . (2) 

This Lagrangian leads to the Heisenberg equations of motion: 

x + {m^ + e^A^)x = (3) 

A + e^x^ A = . (4) 

The Hamiltonian is 

H=lp' + lu\ + l{m^ + e'A^)x\ (5) 

where p{t) = x{t) and H^ = A{t). We next assume that we are in an experimental 
situation where the expectation value of A{t) is so large that quantum fluctuations 
may be ignored. That is, we assume that A is in a classical domain or in a coherent 
state (with large displacement). This is a particular assumption about approximat- 
ing the expectation values involved in taking expectation values of the Heisenberg 
equations of motion, namely: 

{A^x) = {Af {x) ; (x^A) = (x^) (A) (6) 
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Here expectation value means taking a trace with respect to an initial density 
matrix defined at the initial time to, which we take to be to = 0. Taking expectation 
values of the Heisenberg equations with the above factorization, we obtain 

{x) + {m^ + e^ {Af) {x) = 0, (7) 

[A) + e' (x') {A) = 0. (8) 



In this approximation, the equation for (x) is that of an harmonic oscillator with a 
time dependent mass, im?{t) = m? + e^ {A) . A can be thought of as a "classical" 
oscillator (since we do not include its quantum fluctuations) whose mass is determined 
by the quantum fluctuation of the x oscillator (we consider the case where (x)=0). 
The problem of a quantum harmonic oscillator with a time dependent mass can be 
solved in terms of the (numerical) solution of an auxiliary classical oscillator problem. 
We begin by noting that in the Heisenberg picture 

[x{t),p{t)]=t. (9) 

This commutation relation can be satisfled at all times by introducing time- independent 
(deflned at t = 0) creation and destruction operators, a and a\ and using the Ansatz 

x(t) = /(t)a + r(t)at, (10) 

with f[t) satisfying the Wronskian condition 

^[r{t)f{t)-r{t)f{t)] = l. (11) 

The destruction and creation operators a and a) satisfy the usual commutation rela- 
tion [a, a^] = 1. The commutation relation (P) then follows automatically. 

It is easy to show using (H) and (|TU|) that f{t) satisfles the equation of motion 



f + {m' + e'A')f = , (12) 



with the normalization flxed by the Wronskian condition (|TT|). We can either solve 
this classical equation directly numerically, imposing the Wronskian condition at time 
t = 0, or we can automatically impose the Wronskian condition by the substitution 



f{t) = exp 

'2n{t) 



i / Q{t')dt' 







(13) 



where Q{t) satisfies the nonlinear differential equation 
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with 

uj''{t)=m^ + e''A^{t) . (15) 

For simplicity we choose the initial state vector at t = to be the ground state of the 
operator h = a^a, i.e., |^l/(0)) = |0), where a|0) = 0. Then, from ([ToD , the average 
(classical) value of x(t) and p{t) is zero for all time, i.e., {x(t)) = and {pit)) = 0. 
This initial condition pertains in certain semiclassical time evolution problems, such as 
particle production by strong electromagnetic or gravitational fields. For the electric 
field problem, A corresponds to the electric field, and x to the A; = mode of the 
charged particle field (See, e.g., Refs. |]8|l [pJ]] ). 

The quantum fluctuations of x{t) are non-zero and are given by the variable G{t), 



G{t) = {x\t)) = \mf = -^ (16) 



2fi(t) ■ 



From (11^), it is easy to show that G{t) satisfies 



In addition, we find that 

*^W>4(f4). (18) 

The expectation value of Eq. (^ becomes a new effective Hamiltonian 

Heff = {Hit)) 

= lul + 2UlG+^ + lim' + e'A')G. (19) 

The momenta conjugate to G and A are 

Hg = ^ , IIa = A. (20) 

This classical Hamiltonian determines the variables, G and G, necessary for a com- 
plete quantum-mechanical description of the x oscillator. Hamilton's equations then 
yield 

fl^ = -e^AG , (21) 



or equivalently. 



1 fG\ 1 (gV 1 



2 \G) 4 \G) 4G2 
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A + e^GA = 0, 

+ uj^ = , (22) 



which correspond to (|T^ and the expectation value of Eq. (H). 

The classical effective Lagrangian corresponding to the effective Hamiltonian ( |T9|) 
is 

This Lagrangian could also have been obtained using Dirac's action, 

r = j dt{m{t)\i^^ - H\m{t)) ^ jdtL^ff, (24) 

and a time-dependent Gaussian trial wave function as described in Ref . [Q . This vari- 
ational method has recently been used to study the quantum Henon-Heiles problem 
in a mean-field approximation 0. In this method, the Gaussian trial wave function 
is parametrized as follows 

vl>(t) = [2iTG{t)]-'"eM-{^-q{t))\G-\t)/A-iYia{t))+ip{t){x-qm- (25) 

Here, G{t) and T^cit) are the time dependent real and imaginary parts of the width 
of the wave function. One can prove for our problem that if the quantum oscillator 
wave function starts at t = as a Gaussian, it is described at all times by the above 
expression, where G{t) and HgI^) are totally determined by solving the effective 
Hamiltonian dynamics. (For our special initial conditions p{t) = q{t) = 0). Thus 
we find that our effective Hamiltonian totally determines the time evolution of the 
quantum oscillator. 

One interesting "classical" variable is the expectation value of the time dependent 
adiabatic number operator, which corresponds to the number of quanta in a situation 
where the classical A oscillator is changing slowly (adiabatically). For the related field 
theory problem of pair production of charged pairs by strong electric fields (where 
A corresponds to the classical electromagnetic field and x to the k = mode of the 
charged scalar field) this corresponds to the time dependent single particle distribu- 
tion function of charged mesons. To find the expression for the number of quanta, 
which requires the definition of an adiabatic number operator, we begin with the 
wave function for the quantum oscillators corresponding to a slowly varying classical 

background A: 

I r /-t . 1 

(26) 



1 
g{t) = — — exp 

'2uj{t) 



-i / uj(t')dt' 
Jo 



in terms of which we can decompose the quantum operator via 

x{t)=g{t)b{t)+g*{tp{t). (27) 

Requiring the momentum operator to have the form 

p{t)=x{t) = g{t)m+9*{tp{t) (28) 



by imposing g{t)h{t) + g*{t)W(t) = 0, and recognizing that g{t) and g*(t) satisfy 
the Wronskian condition by construction, one finds that b{t) and 6^(t) have the 
usual interpretation as creation and annihilation operators, i.e., [x{t),p{t)] = i and 
[b{t),b^{t)] = 1. Note also that 

b{t)=t[g*{t)x{t)-g*{t)x{t)]. (29) 

It turns out that b^{t)b{t) can be interpreted as a time- dependent number operator 
(assuming a slowly varying (adiabatic) classical field A). The time independent basis 
and the time dependent basis are both complete sets and are related by a unitary 
Bogoliubov transformation, b{t) = a{t)a + /3{t)a^ , where 

a{t) = t[g*{t)f{t)-g*{t)f{t)], (30) 

m = t[g*it)rit)-g*it)rit)], (31) 

and where |a(t)p— |/3(t)P = 1. If we choose as initial conditions, fl{0) = a;(0), Cl{0) = 
ci;(0), then a(0) = 1 and j3{0) = 0. These are the initial conditions appropriate to 
the field theory problem of pair production. The average value of the time-dependent 
occupation number is given by 



n{t) = (b\t)bit)) = \m\' = m^) 



(-")^-3(s-!' 



(32) 



Eq. (|32D allows us to compute the average occupation number of the system as a 
function of time. 

2 Numerical Solution of Hamilton's equations 

We now summarize some previous results from the numerical solution of the Hamil- 
tonian equations obtained in the Gaussian approximation Pj. For calculational pur- 
poses, it turns out to be convenient to scale out the mass via the transformations 
t -^ m^^t, A -^ m~^/'^A, G -^ m'^^G, and e — > ew?!"^ . Then the scaled equations of 
motion are 



1 (G\ 1 (G\^ 1 



2\G A\G 4G^ 



A + e^GA = 

+ 1 + 6^2 = . (33) 



In order to explore the degree of chaos as a function of (scaled) energy and coupling 
parameter e, we calculated surfaces of section and Lyapunov exponents. The surface 
of section is a slice through the three-dimensional energy shell |jl[]. That is, for a fixed 
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Figure 1: A plot of the surface of section for energy = 0.8, e = 1.0, and A = 0. The one 
chaotic region is in the center of the plot. 



energy and coupling parameter, the points on the surface of section are generated as 
the trajectory pierces a fixed place {e.g., A = 0) in a fixed direction. The hallmark of 
regular motion is the cross section of a KAM torus which is seen as a closed curve in 
the surface of section. The hallmark of chaotic motion is the lack of any such pattern 
in the surface of section. In Fig. 1 we show a plot of a surface of section at ii^ = 0.8 
and e = 1 where regular and chaotic regions co-exist. 

The Lyapunov exponent provides a more quantitative, objective measure of the 
degree of chaos. The Lyapunov exponent. A, gives the rate of exponential divergence of 
infinitesimally close trajectories 0. Although there are as many Lyapunov exponents 
as degrees of freedom, it is common to simply give the largest of these. For regular 
trajectories A = 0; for chaotic trajectories the exponent is positive. To define the 
notion of a Lyapunov exponent one begins by considering the infinitesimal deviation 
from a fiducial trajectory. 



Z(t) = lim z(zo + 6,t) — z(zo, t) 
5— >o 



(34) 



where z(zo,t) is a point in phase space at time t with initial position zq. The time 
evolution for Z(t) is given by 



Z(t) = VF 



z(zo,t) ■■^It; , 



m 



where 



z(t) = F(z(t),t) 
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(35) 
(36) 



are the full equations of motion for the system. The (largest) Lyapunov exponent is 
defined to be 



A = lim - In 



m 



z(o) ■ '^^) 

Appendix A of Ref. [^ provides an explicit algorithm for the calculation of all the 
Lyapunov exponents. Since we cannot carry out the t ^ oo limit computationally, 
the regular trajectories are those for which A(t) decreases as 1/t, while the chaotic 
trajectories give rise to A(t) that is roughly constant in time, as judged by a linear 
least-squares fit of log[A(t)] vs. log(t). 

We computed the Lyapunov exponents for three values of the scaled coupling 
constant e (0.1, 1.0, 10.0) and for energies from 0.5 to 2.0. E = 0.5 is the lowest 
energy possible, corresponding to the zero point energy of the oscillator; there is no 
upper limit on E. Fifty initial conditions were chosen at random for each energy bin 
of width 0.1 and coupling parameter. One relevant quantity to study is the chaotic 
volume, the fraction of initial conditions with positive definite Lyapunov exponents 
(corresponding to chaotic behavior). Errors in this quantity arise because of the 
finite number of initial conditions chosen, and because the distinction between zero 
and positive exponents cannot be made with certainty at finite times. We found that 
for e = 0.1, more than 95% of trajectories were regular for all energies tested; for 
e = 1.0 and 10.0, there is a steadily increasing fraction of chaotic orbits between 
0.5 < E < 1.25. For 1.25 < E < 2.0, more than 90% of these orbits are chaotic. In 
Fig. 2 we show that the occupation number is sensitive to initial conditions. 

3 Quantum Oscillators and the Large-N Expan- 
sion 

3.1 The Large-N Expansion and Semiquantum Chaos 

In this section we show that the previous system of equations for the equations of 
motion for two oscillators are just the first term in a large- A^ expansion for expectation 
values of the operator equations of motion of a quantum system consisting of A^ copies 
of the original x oscillator [x ^ Xi,i = 1,2, ...,N) and a single quantum dynamical 
variable A. Such a system is described by the operator equations of motion: 

Xi + {m^ + e^A^{t))xi = 0, 
A + e'(j:xm^A = 0. (38) 




time 



Figure 2: A plot of the occupation number given by Eq. (^) for energy = 1.8, e = 1.0, 
A{0) = 0, Ug{0) = 0. The sohd line is for for G(0) = 0.5; the dashed line is for G(0) = 
0.5001. This plot shows the sensitivity to initial conditions. 



If all the N quantum oscillators Xi have the same initial conditions then, at the level 
of expectation values, we can set all of them equal {{xi) = (x)) and obtain equations 
of motion for the expectation values: 



(x) + m'^ (x) +e'^(^A'^{t)x 
[A) + e^N(x\t)A 



0, 
. 



(39) 



We are interested in initial conditions where (x) = 0. In the large N limit we will 
show that in this case 



x^A 



x^) {A) + 0{l/N) ■ (a^x) = {Af (x) + 0{1/N), 



so that after the rescaling. 



A 



NA 



e'/N 



(40) 



(41) 



we recover the equations for the expectation values (x^) = G{t) and (A) that per- 
tained to semiquantum chaos. To compare with the orginal system describing semi- 
quantum chaos we must solve the quantum system at different N for fixed e. This 
value of e must then be set equal to the value of e used in the semiquantum chaos 
problem. We must also compare (A{t)) to the classical oscillator motion A{t) to see 
the effects of the quantum fluctuations. 
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The exact quantum problem for the coupled quantum oscillators x and A can be 
studied as a function of A^ once we supply the initial wave function at time t = 0. 
Since we are interested in a comparison with our previous calculation, where e is kept 
fixed, it is convenient to write the equations of motion as 

'A\t) 



(x) + m (x) 



g2 



-X 



N 
e-{x^{t)A 



0, 
. 



'A)+eUx'(t)A) = 0. (42) 

The initial conditions, as we change N, are that for fixed Aq{0) at time t = 0, 

{Am = v^io(O) ; (x(0)) = . (43) 



In the Schrodinger picture we have instead the time dependent Schrodinger equa- 



tion 



.d'^{xi,A,t) 
' di 



d 



1 d 



,E^-,,^ + (-^ + e^^^)E 



dx^ 2 dA 



%h A 



^(x„A,t) 



(44) 



With an eye to convenience in solving the Schrodinger equation we will choose an 
initial wave function which is a Gaussian in A and x, so that the initial expectation 
values of A and x will be the ones specified above. This can be implemented by the 
initial wave function: 

N 

*(t=o) = n^x,(^=o)^A(t=o), 



j=i 



^^(0) = [27rD(0)]-i/^exp[-(A-v^io)2(D-i(0)/4-2nz>(0))+2pA(0)(A- 
with 



D(0) + NAl = (A 



t=o 



IId{0) 



m 

AD{0) ' 



Pa{0) 



—I 



-) 



and 



^,,(0) = [27rG'(0)]-i/^exp [-x^(G'-i(0)/4 - ^0^(0)) 



(45) 

NAo)] . 
(46) 

(47) 
(48) 



These initial conditions can be used to compare, as a function of A^, the exact quantum 
problem with the semi quantum chaos problem as well as the 1/N correction to the 
semiquantum chaos problem. 

The large-A^ expansion is best formulated using path integral methods for the 
generating function of the expectation values. For initial value problems rather special 
boundary conditions must be placed on the Green's functions to insure causality. 
The formalism for doing this is the the closed time path (CTP) formulation of the 
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effective action 0. The marrying of the large N expansion to the CTP formahsni 
was accomphshed recently, as described in Ref. |Tl . 

Let us first ignore the issue of boundary conditions on the Green's functions and 
discuss the generating functional for the expectation values. To obtain the generating 
functional we add sources to the original Lagrangian and consider the action 

S = f Ldt (49) 

J c 

where the Lagrangian is, 

L = l^" + T. ir" - l(^" + e'A^)x]+j,x^ + J A . (50) 

The contour c will be chosen in a way that enforces the correct boundary conditions 
for taking expectation values of operators at an initial time t = 0. This will be 
discussed in the appendix. 

The generating functional for the expectation values is given by the path integral: 

Z[J, j] = U[A] /d[xi] . . . U[xj,] e^^[^'-^^'^1 . (51) 



Since (|50|) is quadratic in the Xj variables, we may integrate over all of them in (^ 
and obtain an effective action, given by: 



Z[J,j] = /d[A]e'^^^^[^'^'^'] (52) 

where 

5e//[A;J,j] = l^dt)^--A^ + JAy—TTHG-\A)] 

+ lldt Jjt'J2{Mt)G{t,t%{t')} . (53) 

Here, we have defined 

Go\t,t'; A) = 1^ + (m^ + e'A')\ 6{t - t') . (54) 

Thus the Green's function GQ{t,t') obeys 

1^ + (m^ + e'A')\ Goit,t") = Kt - t") . (55) 

The boundary conditions on this Green's function needed to insure causality will be 
discussed later. We now consider the particular situation with all the Xi identical 

11 



{i.e., have identical initial conditions) so that ji = 
and e as follows: 

A = A/Vn 
J = J/^ 

e = eVN . 



J: -^i 



X. We also rescale A, J 



(56) 

The effective interaction now becomes proportional to A^ as long as e^ is kept fixed 
when A^ is changed. This allows the evaluation of the remaining path integral over 
A by the method of steepest descent and leads to an expansion of the expectation 
values as a power series in 1/A^. The value of the path integral at the stationary 
point is the leading term in this expansion, and the Gaussian fluctuations about the 
stationary point give the 1/A^ correction. Expanding the effective action about the 
stationary point: 



Seff[A;J,j] 



Seff[Ao;J,3] 
+ I dt' ^^^" 



+ 

+ ■ 



dt' 



5A{t' 
dt" 



{A{t') - A,{t')) 



Ao 



6'S, 



ff 



6A{t')6A{t") 



(Aif) 



Ao 



The field Aq is determined by the requirement 



Ao 



and setting, 

the path integral in Eq. 

Z[J,J] 
W[J,j] 



6^S, 



eff 



JAo 



_6A{t) 

^o\t^t') = - sA{t)6A{t') 
I), including terms up to 1/A^, is given by 

Seff[Ao; J,j] + '-TT\n[D^\xo,Ao)] + ■■■ . 



Aoit')) (Ait") - Aoit")) 

(57) 

(58) 



(59) 



(60) 



Since the first term in the action is proportional to A^ and the Gaussian fluctuation 
contribution is of order iV°, the fluctuation term gives the 1/A^ corrections. An 
auxiliary quantity which allows the direct determination of one particle irreducible 
vertices such as inverse Green's functions is the effective action functional F (not to 
be confused with Sejf) which is a Legendre transform of iy[J, j]. Changing variables 
from J,j to the expectation values (x), (A) where 

SW , .. . 6W 



(x) = X 



Sj 



{A)^A 



6J 



(61) 
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We now define the effective action functional (omitting the overhne for symplicity of 
notation) as 

r[a;„ A] = W[j, J]- j dt [mx.,{t) + J{t)A{t)] , (62) 

Jc 

which turns out to be, at order 1/A^, 

■i M i 

r[x„ A] = S,i[x,, A] + — Tr HG,\A)] + -Tr HD,\x, A)] (63) 

where S^ is the classical action and where A and x are the expectation values of the 
Heisenberg operators accurate to order 1/N . 

We now demonstrate that at the stationary phase point the original problem is 
recovered. When the sources are set to zero, 



5S, 



eff 



_6A{t) 
or after rescaling: 



JAo 



*5 + ' 



N 



Y.4iit) + -Go{t,t;A 



Aoit) = , (64) 



dt2 



+ i' 



xl{t) + -Go{t,t;Ao) 



Mt) = 



where Go(t,t'; Aq) satisfies 



— + (m^ + e'A',) Go{t, t'; Ao) = 6{t - t') 



and XQi{t) satisfies 



dt2 



With the identification 



+ (m^ + eM^) Xoi(t)=0 . 



-Go{t,t) = G{t) 



(65) 



(66) 



(67) 



(6^ 



where G{t) = {x'^{t)) — {x{t)) , we arrive at the equations we studied earlier pertaining 
to semiquantum chaos. (For that problem we chose {x{t)) = 0.) Thus we have shown 
that the lowest order in 1/A^ solution to the problem displays semiquantum chaos 
(note that the quantity e^A"^ is invariant under our rescaling). 

In the above, the inverse propagator for the A variable is given by: 



D^\t,t';Ao 



6^S. 



eff 



6A{t)6A{t') 



dQ\t,t'-Ao) + UQ{t,t'-Ao) 



(69) 



^0 



where 



J2xUt) + -Go{t,t;A 



Uo{t,t';Ao) = 2Ne^Ao{t)noit,t';Ao)Ao{t') 



N 



5{t - 1') 



TToit, t'; Ao) = I Goit, t'; Ao) G'o(t', t; Ao) - -J2 ^MGoit, t'; Ao)xoiit') . (70) 
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Rescaling we have: 



do\t,t';Ao) 

Uo{t,t';Ao) 
Mt,t';Ao) 



dt2 



+ e' 



x\t) + -Go{t,t;Ac 



6{t - 1') , 



2g%(t)7ro(t,t';Ao)Ao(t') , 

I Go{t, t'- Ao)G'o(t', t- Ao) - 2x(t)Go(t, t'; Ao)x{t') 



(71) 



We need to invert (p9D subject to the correct causal boundary conditions to find 

DQ{tl,t2]AQ). 

The causal Green's function for initial value problems can be expressed in two 
ways, either as a two dimensional matrix Green's function, or as a path ordered 
Green's function defined on a complex contour. In this paper we will use the second 
method. As discussed in the appendix, the quantity that takes the place of the usual 
Feynman Green's function is the causal Green's function. We begin by defining an 
initial density matrix at time t = by p, and then introducing the two Wightman 
functions (j> and (?< via 



G>(t,t') = i{{x{t)x{t'))-{x{t)){x{t'))} 
G^{t,t') = i{{x{f)x{t)) - {x{t')) {xit))} 



(72) 



where {x{t)x{t')) = Tr{px{t)x{t')}. Time integrals are then defined on the contour 
shown in Fig. 3, with the integration path given by 



dt 



OO ("OO 

dt- dt 

0:c+ JO:c„ 



(73) 



The causal Green's functions which embody the correct boundary conditions are 
then 

G{t,t') = Q,{t,t')G^{t,t') + Q,{t',t)G^{t,t') , (74) 



&c{t,t') 



(75) 



where 

' Q(t, t') for t on c+ and t' on c+ 

for t on c+ and t' on c_ 

1 for t on c_ and t' on c+ 
, 0(t', t) for t on c_ and t' on c_ 

The Green's functions are symmetric and y4>(t,t') = y4<(t',t). These propagators 
take the place of the usual Feynman ones for initial value problems. To calculate 
the 1/A^ corrections to the quantities {A{t)) and (x(t)) and (x^(t)) one needs to 
take functional derivatives of the generating functional, all correct to order 1/A^. 

Specifically we have: 

6W[J,j] 



{m) 



5J{t) 



(76) 
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Figure 3: Complex time contour C for the closed time path propagators. 



where iy[J, j] is given in Eq. (|60D. 

To separate out the leading and next to leading order terms, we write, 

(Ait)) = Aoit) + A,{t) , (77) 

where Aoit) is the quantity determined in first order, by Eq. (pB]). Ai{t) is given by 

'6Do\h,t2;Ao) 



Mt) 



dt' [ dti f dt2Do{t,t';Ao 

J c J c 



Mo(t') 



Do(t2,ti;^). (78) 



By straightforward differentiation, it is possible to show that Ai (t) gets contributions 
from three terms when (x) = 0. Namely, 

Ai(t) = 4")(t) + 4')(t) + 4'=)(t) , 



(79) 



where 



• ~4 



K^''Kt,t2,ts) 
A'-^it) 



Do{t, t2] Ao)Ao(t2)7ro(t2, h] Ao)Do{ts, h; Aq) 

dt2 dt3K^'\t,t2,h) , 



(80) 



i) 5 



Doit, t2] Ao)7ro(t2, h- Ao)Do{t2, h; Ao)Ao{ts) , 
-^ J dt2 J dts J dhK^'\t, t2, ts, ^4^ 

Do{t, t2; Ao)Ao{t2)Go{t2, h; Ao)Ao{t^ 
xGo{t2,U; Ao)Ao{U)i:{ts,U; Ao) . 
15 



^2) 



In the above expressions, the self energy S is given by 



S(t3, ti] Ao) = iDo{t3, U] v4o)G'o(^3, ^4; ^0 



and the polarization vr by 



7r(t3, t^] Aq) = iGo(t3, ti, Ao)Go{t3, t^, Aq 



^3) 



U) 



Note that the causal Green's functions are symmetric under the interchange of time 
labels. It is easy to show, using the rules for convoluting the causal Green's functions 
(given below), that Ai{t) only depends on information from previous times. 

In order to evaluate the above graphs we first have to determine the causal prop- 
agator Dq. From (B9), we have 



D,\t,t')=d^\t,t') + Uo{t,t') 



^5) 



where cig (t, t') and no(t, t') are defined in ([701). The defining equations for the inverses 
are 



If we now put 
then (RBI) implies 



dt'do\t,t')dQ{t',t") = 6{t-t") , 
j dt' D^\t,t')D^{t',t") = 6{t-t"). 

J c 

Doit,t') = doit,t') + Doit,t') , 
Jdt' {do\t,t') + no(t,t')} ^o(^'>^") = ^(t - 1") 



i7) 



or 



J^dt'do\t,t')D,{t,t') = -J^dt'U,{t,t')D,{t',t") 
We can now invert the left hand side using (|86D again, obtaining: 



(89) 



Do{h,U) = doituU) - Jjt2 Jjt3do{tuh)Uo{t2,t3)Do{h,U) . (90) 

We use (pO|) to find the inverse. Note that do{t,t') satisfies 



dt2 



+ e^ 



A^ 



Y^xUt) + -Go{t,t) 



do{t,t') = 6{t-t') . 



(91) 



We will need to solve for this Green's function c/q in order to find Doit^t') from 
(pop. As shown below this can be done in analogy to the determination of G{t) in 
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the semiquantum problem, that is, we will introduce a set of mode functions for an 
auxiliary quantum problem. 

To study the time it takes for the 1/A^ corrections to become significant, we also 
need to determine the 1/A^ correction to the x oscillator propagator. This is obtained 
once again by functional differentiation. Now we turn to the full Green's function 
G{t,t') for the x oscillator (to order 1/N), from which we can obtain 

(^x\t)) = -G{t,t). (92) 

This is most easily determined to order 1/A^ from the inverse Green's function which 
by the chain rule is the negative of the second derivative of the effective action func- 
tional r[a;, A] with respect to x. The effective action functional F to order 1/A^ is 

AT i 

T[x,,A] = S4x,,A\+t-TT\n[G,\A)] + -TT\n[D^\x, A)] , (93) 

where Sd is the classical action and here A and x are the expectation values of the 
Heisenberg operators accurate to order 1/A^, i.e., the A in this equation is the solution 
of 6T/SA = 0, or 






J2xl,{t) + -Go{t,t;A) 

t 



iff I 

A{t) = -J^dt^J^dt2Do{ti,t2;A)- 



6A{t) 
(94) 

and differs from the stationary point Aq of the original path integral by terms of order 

1/N. Now 

^ f ,. f ,. r^ u . ,J^Doit2,h;A) 



Go\t,t';A)6ij-'-J dtij dt2Do{ti,t2;A)- 



6xi{t)5xj{t') 
6ijG-\t,t';A) . (95) 



And, since 

2e'6{ti-t2)6{ti-t)6{t2-t') 



6^Do{t2,h;A) _ ^ 2, 



6xi{t)6xj{t') 

- 4e^A{t)Go{t,t'; A)A{t')[6{h - t)6{t2 - t') + 5{t^ - t')6{t2 - t)] , 

(96) 

we obtain 

G'\t, t'- A) = Go-\t, t'- A) + Gf\t, t; A), (97) 

where Gi is of order 1/N and given by two terms. Explicitly, 

~2 ~4 

G,-\t, f) = -i^D.it, t)5{t - t') + 4^io(t)S(t, t')A,{t') . (98) 
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To obtain the 1/A^ expansion of G we have to invert this to order 1/A^. We first need 
to reexpand Gq^{A) up to order 1/A^ since A = Aq + Ai where Ai is given by ([78|) . 
We have, 

G-\t,t') = G'o"'(Ao)(t,t') + 2e%(t)Ai(t)(5(t-t')+Gr'(t,t') 

^ Go-^(Ao)(t,t') + ^A(t,t'). (99) 

Inverting to order 1/A^ we obtain finally 

G{t, t') = Go{t, t'- Ao)-^l dti ^ dh Go{t, ti; Ao)A{h,h)Go{t2, t'; Ao) + 0{1/N^) . 

(100) 

3.2 Closed time path contour and causality 

In this section we discuss the causality of the Green's functions as given by the CTP 
formalism. We will evaluate the time integrals using the closed time path (CTP) 
contour shown in Fig. 3. The integration path is given explicitly by 

/• /"OO /"OO 

dt= dt- dt . (101) 

Jc JO:c+ JO:c_ 

The Green's functions are now given by functions of the form, 

A{t, t') = 0,(t, t')A> {t, f) + ee(t', t) A< (t, t') , (102) 



where 0c(t,t') is defined in Eq. ([751). We call such functions causal. The causal 
Green's functions are symmetric implying thereby that A>(t, t') = A^{t',t). In order 
to prove causality of any particular graph we need to discuss two lemmas. The first 
is that if we have a loop of two causal functions, such as found in a self energy graph, 
then that is also a causal function. To show this one just needs the definition of Q^- 
If the two causal functions are 

B{t,t') = Q,{t,t')B^{t,t') + Q,{t',t)B4t,t') , 

G{t,t') = Q,{t,t')CUt,t') + Q,{t',t)C<{t,t') , (103) 

then the self energy graph is also causal. Letting the self energy 

A{t,t') = iB{t,t')C{t,t') , (104) 



and setting. 



we find 



A>,<(t,t') = ^5>,<(t,t')C>,<(t,t') , (105) 

A{t, f) = e,(t, t')Ay{t, f) + 0e(t', t)A<(t, t') , (106) 
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which is the desired result. 

The next lemma is that the matrix product of two causal functions is causal. That 
is, if we think of the Green's functions as matrices in time and if we matrix multiply 
B and C so that A{ti,t3) is given by 

A{t,,t3)= J^dt2B{tut2)Cit2,h) , (107) 

we find then 

A{t,t')=e,{t,t')A^{t,t') + &,{t',t)A^{t,t') , (108) 

where 

A>{tuh) = - f'dt2B>{ti,t2)[C>{t2,h)-C^{t2,h)] 
< Jo < 

+ r dt2[B^ih,t2) - B^{t,,t2)]C>it2,h) . (109) 

JO < 

This lemma is discussed in Ref. |Tl| and is obtained directly by breaking the time 
integration into three segments, viz., 

(i) <ti <t, on C+ , 

(a) t <ti < oo, on C+ , 

(iii) <ti < oo, on C_ . (110) 

One then uses the definition of the G function (|75D and collects all the non-cancelling 
terms. 

Now consider the product of three causal functions: 

Aih,U) = I dt2 I dhB{ti, t2)C{t2, h)D{h, U) . (Ill) 

J c J c 

We can work this case out by applying the second lemma from left to right. That is, 
we can let 

E{t^,h) = j^dt2B{h,t2)C{t2,h) . (112) 

Then E{ti,t-i) is causal and is given by an equation of the form ( |108| ). We are then 
left with an equation of the form: 

A{ti,U) = f dUE{ti,h)D{h,U) , (113) 

J c 

which is also causal. In the same way, we can find causal relations for any number of 
CTP integrals of causal functions. We can also apply the lemma from right to left. 
After doing the integrals sequentially one is eventually left with 

f{t)= fdhF{t,h)= /*[F>(t,ti)-F<(t,ti)] , (114) 

Jc Jo 

which exphcitly displays the causality. 
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3.3 Lowest order causal Green's functions 



With the resuhs of the previous section in hand we are now in a position to solve for 
the Green's functions with the correct causal structure. We would like to solve 



— + K + e2i^(t)) Go(t,t') 



d^ 
dt2 



e^ 



xl{t) + -Go{t,t) 



do{t,t') 



S{t - 1') , 
6{t - 1') , 



;il5) 



subject to causal boundary conditions and an initial density matrix which is that of an 
adiabatic vacuum at time zero. The adiabatic requirement is satisfied by considering 
auxiliary quantum fluctuation operators Xq and A^g such that {xq) = (^og) = 0. (Note 
that we are defining the fluctuation operators by writing the Heisenberg operator 
A = {A) + Aq.) These operators are defined via 



xoit) = f{t)a + r{t)a* 
Ao,{t) = g{t)h + g*{t)h* , 

where a and h are canonical annihilation (creation) operators satisfying 

[a, a^] = [b, b^] = l , 

and the adiabatic vacuum is defined by 



a|0)., = &|0) 



ad 



. 



;ii6) 



;ii7) 



:ii8) 



The /(t) and g{t) are functions of time satisfying the homogeneous equations, 



^, + {m' + e'Al{t))\fit) 



dt2 



+ e' 



xlit) + -Go{t,t) 



git) 








;ii9) 



with Wronskian conditions, 



4r(t)A/(,)} 



t\g*{t)-git) 



1 . 



(120) 



We can write the causal Green's functions in terms of the complex functions f{t) and 
g(t),as follows: 



Go{t,t') = e,{t,t')Go>it,t') + e,it',t)Go4t,t') 

doit,t') = 0e(t,t')c?>(t,t') + ec(t',t)rf<(t,t') , 
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(121) 



where 



Go>(t,t') = (xo(t)xo(t'))=^/Wr(t'), (122) 

Go<it,t') = {xoit')xoit)) = tfit')rit), (123) 

d>(t,t') = {Ao,{t)Ao,{t'))=zg{t)g*{t') , (124) 

d<(t,t') = {Ao,{t')Ao,{t))=tg{t')g*{t). (125) 



Note also that 



Go{t,t) = t\f{t)\' = tG{t), 
do{t,t) = t\g{t)\\ (126) 

Here the factors of i are introduced to agree with earher definitions. We are now in a 
position of being able to find the 1/A^ quantum corrections to Ai{t) and {x'^{t)), using 
Eq. (pOD to construct DQ{t,t'). The last-named quantity obeys the integral equation 

Do{ti,U) = doituh) - fdt2 f dtsdo{tut2)Ilo{t2,h)Do{h,U) • (127) 

Jc J c 

We want to evaluate this in terms of the lowest order quantities. The polarization 
can be put in causal form since 

7ro(t, t') = i[Q,{t, t')Gl^{t, t') + e,(t', t)Gl^{t, t')] (128) 

and 

no(t,t') = 2g^^(t)7^o(t,t')^(i') • 

= ee(t,t')n>(t,t') + 0c(t',t)n<(t,t') (129) 



Now we can do the integrals sequentially using ( p.09| ). First we do the integral 
over ^2- Then we write Dq schematically as 

Do{t,t') = do{t,t')- f dtiUit,h)Do{h,t') 

Jc 

= do{t,t') - D{t,t') (130) 

where U{t,ti) is causal and is determined by the matrix multiplication of do and 11. 
Then 

D>{t,t') = - f dt^U>{t,t^)[Doyih,t')-Do^ih,t')] 
< Jo < 

+ f dh[U^{t,t,) -U^{t,t,)]D >m,t') , (131) 

JO o< 
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with boundary conditions at t = t' = 0: 

Do(0,0) = rfo(0,0); D{0,0) = 0. (132) 

In order to determine the causal matrix (in time) D{t,t') one recognizes that for 
causal Green's functions 

Ay{t,t') = A<(t',t), A^{t,t') = -A^{t',ty . (133) 

When doing numerical work, the time integrals are replaced by discrete sums, 
with t = me and t' = ne where e is the time step. The explicitly causal update is 
then, for m > n, 

m— 1 

Doy{m,n) = ^[t/>(m, /c) -f/<(m, A;)]L'o>(^,^) 

fc=0 

n-l 

- Y.lU>im,k)[DoUk,n)-Do^ik,n)]. (134) 

fc=0 

Thus starting with the known value of Dq{0,0) we can construct the entire causal 
propagator Do(t,t') using this and the relations (|133|) . 



3.4 Initial conditions 

Particle production in the early universe and particle production by a classical electric 
field are two external field problems which admit a particular type of intial condition. 
One starts these initial value problems off with no particles in the appropriate matter 
field and with the external field in a "classical" state (one where the expectation value 
dominates quantum fluctuations). This is also the situation here with A corresponding 
to the external field and the Xi to the modes of the quantum matter field. 

In line with the above desires for the appropriate initial condition, we want to 
enforce {xoi(t)) = {Ag(t)) = for all i and values of t. This is accomplished by 
choosing the initial state to be an adiabatic vacuum a |0)^^ = b |0)^^ = 0. 

Next, we note that solutions of ( |119| ), 

rd_2 ^ 

( ^2 



^^ + {m' + i'Al{t))\f{t) = 



^^ + eVWP|^(t) = 0, (135) 

which automatically obey the Wronskian condition can be written in the form: 

fit) = ^ e-'/o"^^*')'^*' , 

git) = ^ ^-ij;;n,it')dt' ^^3g^ 
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where ^f{t) and Qg(t) are solutions of the nonhnear equations, 

2 






2\n'J^\n'ri +^^ " '^^' 



1^^.^ , ^^V + fiJ ^ ^2^ (137) 



2 V^s/ V^.. 
with 

We want to match our solutions to asymptotic (adiabatic) Heisenberg operators. 
This will be accomplished by making the choices: 

nj{0) = ujf{0); nf{0)=uj{0); fi,(0) = 07,(0); (1,(0) = ^M ■ (139) 

Finally, the following initial conditions for f{t) and g{t) are obtained 

/(o) = -A=; 9(0) ^ 



v/2^/(0) ' ^2^/(0) 

^(Q) ^ ^/(0)/2^/(0)+^^/(0) . .(o)_ c:;,(0)/2a;,(0)+^a;,(0) ^^^^^ 
/2^/(0) ' ^2^.(0) 



These initial conditions have to be supplemented by the initial values for Aq{0) and 

io(o). 

The solution of the 1/A^ correction to {x'^{t)) will tell us how the time scale for 
breakdown of the 1/A^ expansion depends on A^. (In certain semiclassical problems, it 
was found in Ref. [|T^ that the breakdown time went as logA^.) We will also be able 
to determine how this time scale is related to the Lyapunov time scale. The realm of 
validity of the large N expansion is determined by comparing to the exact quantum 
problem which we will solve as a function of A^. For the quantum problem we need 
to supply the parameters of the initial Gaussian. For the initial data specified above 
(adiabatic initial conditions), 

^(^ = 0) = ^T^ ' ^(^ = 0) = ^T^ ' ^141) 

2u}f{0) 2ujg{0) 
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This information determines the real and imaginary part of the width of the initial 
wave function for the exact calculation. We have started numerical simulations of 
both the 1/A^ expansion and the exact quantum problem. We hope to be able to 
present our findings in the near future. 
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A The Closed Time Path Formahsm 

In scattering theory one is interested in the probabihty that an initial state evolves 
into a particular final state. The boundary conditions for the Green's functions for 
the correlation functions in that situation are the Feynman ones, and these correlation 
functions can be obtained from the conventional path integral formalism which defines 
transition elements between states at one time, t (usually taken to be in the infinite 
past) to states at another time t' (in the distant future). If the class of paths is 
restricted to be the vacuum configuration at both of its endpoints, then the two 
states are the \in) and {out\ vacuum states of scattering theory respectively. The 
generating functional Z[J,j] of for the Green's function of scattering theory is the 
transition matrix element 

Z[J,j]{t,t') = {out,t'\in,t)jj (144) 

in the presence of the external sources J and j. 

By varying with respect to the external sources we obtain matrix elements of the 
Heisenberg field operators between the \in) and {out\ states. For this reason we may 
refer to the conventional formulation of the generating functional Z as the "in-out" 
formalism. The time-ordered Green's functions obtained in this way necessarily obey 
Feynman boundary conditions, and these are the appropriate ones for the calculation 
of transition probabilities and cross sections between the \in) and {out\ states. On 
the other hand the off-diagonal transition matrix elements of the in-out formalism are 
completely inappropriate if what we wish to consider is the time evolution of physical 
observables from a given set of initial conditions. The in-out matrix elements are 
neither real, nor are their equations of motion causal at first order in 1/N, where 
direct self interactions between the fields appear for the first time. What we require 
is a generating functional for diagonal matrix elements of field operators with a cor- 
responding modification of the Feynman boundary conditions on Green's functions 
to ensure causal time evolution. This "in-in" formalism was developed more than 
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thirty years ago by Schwinger, Bakshi and Mahanthappa and later by Keldysh, and 
is called the closed time path (CTP) method p. 

The basic idea of the CTP formalism is to take a diagonal matrix element of the 
system at a given time t = and insert a complete set of states into this matrix 
element at a different (later) time t' . In this way one can express the original fixed 
time matrix element as a product of transition matrix elements from to t' and the 
time reversed (complex conjugate) matrix element from t' to 0. Since each term in this 
product is a transition matrix element of the usual or time reversed kind, standard 
path integral representations for each may be introduced. If the same external source 
operates in the forward evolution as the backward one, then the two matrix elements 
are precisely complex conjugates of each other, all dependence on the source drops 
out and nothing has been gained. However, if the forward time evolution takes 
place in the presence of one source J+ but the reversed time evolution takes place 
in the presence of a different source J_, then the resulting functional is precisely the 
generating functional we seek. Indeed (setting j = and A^ = 1 here for simplicity). 



'^in[J+-, 'J-\ 



[Vm]{in\i 
\D'^]{in\T* exp 

(^,t'|Texp 



t' 



-I / dtJ^{t)<i){t) 



t' 



|^,t') X 



I / dtJ+{t)(j){t) 





\in) 



so that, for example. 



6Mt) 



5iy,„[j+,j_ 



J.=J_0 



6J4t) 



J+=J-=0 



{in\(f){t)\in) 



(145) 



(146) 



is a true expectation value in the given time-independent Heisenberg state \in). Here 
(j)(t) = {A(t),x{t)}. Since the time ordering in Eq. (|145|) is forward (denoted by T) 
along the time path from to t' in the second transition matrix element, but backward 
(denoted by T*) along the path from t' to in the first matrix element, this generating 
functional receives the name of the closed time path generating functional. If we 
deform the backward and forward directed segments of the path slightly in opposite 
directions in the complex t plane, the symbol 7c may be introduced for path ordering 
along the full closed time contour, C, depicted in Fig. 3. This deformation of the path 
corresponds precisely to opposite ie prescriptions along the forward and backward 
directed segments, which we shall denote by C± respectively in the following. 

The doubling of sources, fields and integration contours in the CTP formalism 
may seem artificial, but in fact it appears naturally as soon as one discusses the time 
evolution not of states in Hilbert space but of density matrices. Then it is clear that 
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whereas | ) ket states evolve with Hamiltonian H, the conjugate ( | bra states evolve 
with —H, and the evolution of the density matrix requires both. Hence a doubling 
of all sources and fields in the functional integral representation of its time evolution 



kernel is necessary. Indeed, it is easy to generalize the functional in ( |145| ) to the case 
of an arbitrary initial density matrix p, by defining 



Z [J+, J_,p\ = Tr <^ p r* exp 



i / dtJ^{t)(t){t) 



Texp 



i / dtJ+{t)(p{t) 



t' 



t' 
-I I dtJ^{t)(t){t) 
/o 



x{il!\Texp 



I I dtJ+{t)<p{t) 





1^) 



(147) 



Variations of this generating function will yield Green's functions in the state specified 
by the initial density matrix, i. e. expressions of the form. 



Tr{p0(ti)0(t2)0(t3)...} . 



(148) 



Introducing the path integral representation for each transition matrix element in 
Eq. (|147D results in the expression. 



exp 



oo 

II dt{ L[ 





-L\ 



■J+<P+ - J- 



(149) 



where L is the classical Lagrangian functional, and we have taken the arbitrary future 
time at which the time path closes t' —>■ oo. 

The double path integral over the fields 0+ and 0_ in ( |149D suggests that we 
introduce a two component contravariant vector of field variables by 



a = 1,2 



(150) 



with a corresponding two component source vector, 

J'^= I ^+ I ; a = 1,2. 



(151) 



Because of the minus signs in the exponent of ( |149|) , it is necessary to raise and lower 
indices in this vector space with a 2 x 2 matrix with indefinite signature, namely 



Cab = diag (+1,-1) = c"* 
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(152) 



so that, for example 

rcab^' = J+(t)+ - J-<P^ . (153) 

These definitions imply that the correlation functions of the theory will exhibit a 
matrix structure in the 2x2 space. For instance, the matrix of connected two point 
functions in the CTP space is 

5^W 



G''\t,t') 



j=o 



6Ja{t)6Mt') 

Explicitly, the components of this 2x2 matrix are 

G^\t,t') = G>(t,t')= *Tr{p$(t)$(t')}con, 
G'\t,t') = G'<(t,t') = ^Tr{p$(t')$(t)}co„, 

G'\t,t') = tTi{prmt)^t')]}^^^ = e{t,t')Gyit,t') + e{t',t)G^it,t') , 

G^\t,t') = «Tr{p r*[$(t)$(t')]}^^^ = e{t',t)G>{t,t') +e{t,t')G<{t,t') 
Notice that 



(154) 



with the usual convention that 



G^'{t,t) =G'\t,t) 



0{t,t) = - 



(155) 



(156) 



The 2x2 matrix notation has been discussed extensively in the literature 0. 
However, the development of the CTP formalism is cleaner, both conceptually and 
notationally, by returning to the definition of the generating functional ( p.47|) , and 
using the composition rule for transition amplitudes along the closed time contour in 
the complex plane. Then we may dispense with the 2x2 matrix notation altogether, 
and write simply 



[DiP]{ip'\T*exp 



dtJ.{t)(t){t) 



l^)(^l^exp 



dtJ+{t)^{t) 



\^) 



{ip'\Tcexp 



i / dtJ{t)(t){t) 



W) 



(157) 



so that ( |147|) may be rewritten more concisely in the CTP complex path ordered 
form. 



Zc [J. p] 



Tr|p(^Tcexp 


i f dtJ{t)4){t) 
. Jc 


)} 






j[V^']j[V^^] (^V^' )£ Mexp 


i 1 dt (L[0] + J0) 


.(158) 










(159) 



The advantage of this form is that it is identical in structure to the usual expression 
for the generating functional in the more familiar in-out formalism, with the only 
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difference of path ordering according to the complex time contour C replacing the 
ordinary time ordering prescription along only C+. Hence, all the functional formalism 
of the previous section may be taken over line for line, with only this modification of 
complex path ordering in the time integrations. For example, the propagator function 
becomes 



Git.t') = ecit,t')G>{t,t') + ecit',t)G4t,t') 
= ec{t,t')G^\t,t') + ec{t',t)G^\t,t') 

where Oq is the CTP complex contour ordered theta function defined by 

' 9{t,t') for t,t' both on C+ 



(160) 



Oc{t,t') 



9{t',t) for t, t' both on C_ 
1 for t on C_ , t' on C+ 

for t on C+ , f on C^ 



(161) 



With this definition of G{t, t') on the closed time contour, the Feynman rules are the 
ordinary ones, and matrix indices are not required. In integrating over the second 
half of the contour C_ we have only to remember to multiply by an overall negative 
sign to take account of the opposite direction of integration, according to the rule, 

/■ /"OO /"OO 

dt= dt- dt . (162) 

Jc Jo c+ Jo C- 

A second simplification is possible in the form of the generating functional of 
( p,59| ), if we recognize that it is always possible to express the matrix elements of the 



density matrix as an exponential of a polynomial in the fields |T2 

R + Ra{to)ip''ito) + Rab{to)^\to)^\t 



(/IpI^' ) 



exp 



+ ... 



(163) 



Since any density matrix can be expressed in this form, there is no loss of generality 
involved in expressing p as an exponential. If we add this exponent to that of the 



action in ( 159 ), and integrate over the two endpoints of the closed time path ip^ and 
if^, then the only effect of the non-trivial density matrix p is to introduce source 
terms into the path integral for Zc[Jip] with support only at the endpoints. This 
means that the density matrix can only influence the boundary conditions on the 
path integral at t = 0, where the various coefficient functions Ra-, Rab, e.tc. have the 
simple interpretations of initial conditions on the one-point (mean field), two-point 
(propagator), functions, etc. It is clear that the equations of motion for t 7^ are 
not influenced by the presence of these terms at to = 0. In the special case that 
the initial density matrix describes a thermal state, p/3 = exp{— /?i7} then the trace 
over pp may be represented as an additional functional integration over fields along 
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the purely imaginary contour from t = —i(3 to t = traversed before C_ in Fig. 3. 
In this way the Feynman rules for real time thermal Green's functions are obtained 
p!3| . Since we consider general nonequilibrium initial conditions here we have only the 
general expression for the initial p above and no contour along the negative imaginary 
axis in Fig. 3. 

To summarize, we may take over all the results of the usual scattering theory 
generating functionals, effective actions, and equations of motion provided only that 
we 

1. substitute the CTP path ordered Green's function(s) ( |160| ) for the ordinary 
Feynman propagators in internal lines; 



2. integrate over the full closed time contour, C, according to ( |162| ); and 



3. satisfy the conditions at t = corresponding to the initial density matrix p. 
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